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ABSTRACT 

Context. In nature we observe galaxy aggregations that span a wide range of magnitude gaps between the two first-ranked galaxies of 
a system (Amn). Thus, there are systems with gaps close to zero (e.g., the Coma cluster), and at the other extreme of the distribution, 
the largest gaps are found among the so-called fossil systems. The observed distribution of magnitude gaps is thought to be a conse¬ 
quence of the orbital decay of M* galaxies in massive halos and the associated growth of the central object. As a result, to first order 
the amplitude of this gap is a good statistical proxy for the dynamical age of a system of galaxies. Eossil and non-fossil systems could 
therefore have different galaxy populations that should he reflected in their luminosity functions. 

Aims. In this work we study, for the first time, the dependence of the luminosity function parameters on Ami 2 using data obtained by 
the fossil group origins (TOGO) project. 

Methods. We constructed a hybrid luminosity function for 102 groups and clusters at z < 0.25 using both photometric data from 
the SDSS-DR7 and redshifts from the DR7 and the FOGO surveys. The latter consists of ~1200 new redshifts in 34 fossil system 
candidates. We stacked all the individual luminosity functions, dividing them into bins of Ami 2 , and studied their best-fit Schechter 
parameters. We additionally computed a “relative” luminosity function, expressed as a function of the central galaxy luminosity, 
which boosts our capacity to detect differences—especially at the bright end. 

Results. We find trends as a function of Amyi at both the bright and faint ends of the luminosity function. In particular, at the bright 
end, the larger the magnitude gap, the fainter the characteristic magnitude M". The characteristic luminosity in systems with negligible 
gaps is more than a factor three brighter than in fossil-like ones. Remarkably, we also find differences at the faint end. In this region, 
the larger the gap, the flatter the faint-end slope a. 

Conclusions. The differences found at the bright end support a dissipationless, dynamical friction-driven merging model for the 
growth of the central galaxy in group- and cluster-sized halos. The differences in the faint end cannot be explained by this mecha¬ 
nism. Other processes—such as enhanced tidal disruption due to early infall and/or prevalence of eccentric orbits—may play a role. 
However, a larger sample of systems with Anin > 1.5 is needed to establish the differences at the faint end. 

Key words. Keywords should be given 


1. Introduction 

The exi stence of fossi l galax y groups was proposed for the first 
time bv iPonman et alJ (Il994l) . In that work, it was suggested that 
the isolated elliptical galaxy RX J1340.6-I-4018 was probably 
an evolved compact group of galaxies. They claimed that those 
galaxies that were close to the center of the system could have 
merged in a single elliptical galaxy. This is why they are called 
’’fossil groups” (FGs). The most accepted o bservational defini - 
tion for this kind of object was proposed bv iJones et alJ (l2003h . 
They defined a system of galaxies as a fossil if it presents a mag¬ 


nitude gap of at least two magnitudes between the two bright¬ 
est member galaxies (Ami 2 > 2 ) in the r-band within half of 
its virial radius and if the central galaxy is surrounded by an ex¬ 
tended X-ray halo of Lx.boi > 50 ■ The latter criterion 

was adopted to distinguish large isolated galaxies from group¬ 
sized systems, but it is a lower limit, so it does not exclude the 
existence of ’’fossil clusters” (as proposed by I Cvprian o et al.l 


I 2 OO 6 I : [Mendes de Oliveira et alJl200^IZarattini et aljl2014 ). For 

this reason, we refer to fossil systems, but we prefer to maintain 
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the classical notation of FGs, which is usually accepted in the 
literature, as we did in the other papers of this series. 


Different observational properties of FGs were studied. The 
properties of the hot intracluster component were analyzed us¬ 
ing scaling relations that include some X-ray properties of the 
system. The Lx - Tx relation is generally similar to that of nor - 
mal clusters (see lKhosroshahi et al .11200^ iHarrison et al.ll2012h . 
whereas differences have been found in scaling relations that 
combine both optical and X-ray properties, such as the Lx - Topt, 
Lx - cr^, and Tx - <Ty relations. In fact, some works exploring 
these relations suggest that fossil systems are brighter in the X- 
ray range (or fainter in the o ptical range) than n ormal groups 
and clusters (for exampl e, see iProctor et al.ll201 ih . In contrast, 
iKhosroshahi et alJ (12014l) find that fossils are underluminous in 
the X-ray range. Howeve r, these differences can be attributed to 
obser vational bi ases (see Voevo dkin et al.ll201fll : iHarrison et alJ 
l2012h . Recently. iGirardi et al.l ( 2^4l) have analyzed a sample of 
15 spectroscopically confirmed FGs, finding no significant dif¬ 
ferences in the Lx - Topt relation. Particular attention was paid to 
the homogeneity of the data set and the analysis process in this 
work. 


The brightest group galaxies (BGGs) of FGs are considered 
the most massive galaxies in the Universe, and as such they have 
been s tudied well in the literature. For example, IHarrison et alJ 
(12012li show that both the absolute magnitude of the BGGs 
and the fraction of light contained in them are correlated with 
the magnitude gap, r esults that we have recently confirmed in 
IZarattini et alJ (l2014li . Moreover, thei r luminosity is correlate d 
with the system velocity dispersion (IKhosroshahi et al]|2006h . 
Observations of the BGG s isophotal shape are not conclusive. 
IKhosroshahi et al.l (l2006h show that these objects present disky 
isophotes in the central part, whereas both isolated ellipticals 
and c entral ellipticals i n clust ers s how boxy isophotes. I n con- 
trast, iLa Barbera et al.l (l2009h and iMendez-Abreu et ^ (12012l) 
find no differences in this sense between fossil and non-fossil 
systems. In addition, the size-luminosity relation, the fundamen¬ 
tal plane, and the Faber-Jac kson relation are similar for f ossil and 
non-fossil central galaxies (iMendez-Abreu et'^l2012h . Recent 
studies of the stellar population of BGGs seem to indicate that 
their age, metallicity, and a enhanceme nt are similar to those o f 
central galaxies in non-fossil systems JUa Barbera et al.ll200^ . 
Moreover, the absence of large gradients in the metallicity radial 
profiles rules out the hypothesis of the mon olithic collapse for 
BGGs in FGs dEigenthaler & Zeilinged20T^ . In summary, these 
observational properties indicate that BGGs in fossil and non¬ 
fossil systems show similar properties. The only relevant differ¬ 
ence is the fraction of light enclosed in the BGG. This shows that 
BGGs in fossil groups may have formed via similar (but perhaps 
more efficient) physical mechanisms to non-fossil ones. 


All these observational properties can be explained in terms 
of the formation scenario of FGs. Numerical simulations show 
that the halo of a FG comprises half of its mass at z > 1 
( D’Onghia & Lakj l2004t iD’Onghia et alJ l2005t iDariush et alJ 
1201 nh . Then" it grows via minor mergers alone, accreting 
only one third of the galaxies of regular groups or clusters 
dvon Benda-Beckmann et al.l l2008l) . Moreover, in simulations, 
FGs always show an assembled mass that is, on average, highe r 
than non-fossil systems at any redshift dPariush et al.l l2007h . 
These simulations seem to favor what is considered to be the 
“classical” formation scenario for FGs. They are thought to be 
very old systems that were able to assemble the majority of their 
mass at high z, where the M* galaxy population has been canni¬ 
balized by the BGG. 


Nev ertheless, the formation of the BGG could be a long-term 
process. iDfaz-Gimenez et ^ d2008h claim that the last major 
merger for the BGG occurs at a later time in fossil than in non¬ 
fossil systems. iGozaliasl et alJ d2014l) suggest that the BGGs of 
fossil systems are the result of multiple mergers of M* galaxies 
in the past 5 Gyr. Moreover. Ivon Benda-Beckmann et ^ d2008h 
claim that the fossil phase could be only transitional and that 
the interaction with other groups or clusters could erase the gap 
in magnitude. There is another possible scenario that is com¬ 
pletely different from all of those mentioned above. This is the 
so-called failed group scenar io, which was first proposed by 
iMulchaev & ZabludoftI dl999l) . In this scenario the gap in mag¬ 
nitude is not due to the evolution of the system; rather, it appears 
by chance during the formation of the system itse lf. However, 
recen t simulations matching subhalo abundance dHearin et al.l 
12013h seem to indicate that this scenario is not a good repre¬ 
sentation of reality. 

The short formation time described in the classical scenario 
would give fossil systems enough time to merge all M* galax¬ 
ies to form the massive central galaxy. The M* galaxies are the 
natural candidates for the merging process, since the dynami- 
cal friction—which is responsible for the decay of the orbits 
(IChandrasekhailfT943l) —is higher for more massive satellites. 
Moreover, FGs are supposed to host M* galaxies in more radial 
orbits, and this can boost the efficiency of the mer ging process 
(ISommer-Larsenl2006l:lBoylan-Kolchin et al.l2008l) . This means 
that fossil and non-fossil systems should haye different luminos¬ 
ity functions (LFs). 


The LF giyes the number density of galaxies per lumi¬ 
nosity interyal, and it is a yery powerful tool for studying 
the galaxy population in groups and clust ers. A recent re - 
yiew of the principal results can be found in iJohnst^ (1201 ll) . 
In the case of FGs, because there are so few known sys¬ 
tems, the majority of publications haye analyzed the LFs of 
indiyidual FGs. In particular, each analyzed system seems to 
sho w a peculiar LF that does not accord with the others 
(see KhosroshahieLaLl EOO^ |hfendes_dg_01iyeira_gtaL 200^ 


lAguerri et al.l201 1 : Adami et al.l20121: IKhosroshahi et al.' 20141) . 

Thus, a systematic and homogeneous study is still required. For 
this reason, we present a large study here of a sample of 110 sys¬ 
tems, containing 19 confirmed fossi ls. The criteria in the defini- 
tion of fossils are those reported in IZarattini et al.l (l2014l) . With 
this unique data set, we are able to present the first study of 
the dependence of the LF on the magnitude gap in group- and 
cluster-sized systems within half the R 200 radius. 


This work is part of the FOssil Group Origins (FOGO) 
project. This is a multiwayelength study focused on th e sam¬ 
ple of 34 FG candidates proposed by lSantos et al.l (l2007l) . A de¬ 
tailed oyeryiew o f the FOGO project is presented in the first pa- 
per of the series, Aguerri et al.l (|201 ll hereafter FOGO I). In 
IMendez- Abreu et al.l ( 20121 hereafter FO GO II), we explored 
the properties of the central galaxies; in IGirardi et al.l (12014 
hereafter FOGO III), we presented the study of the Lx-Lopt 
relations; and the ch aracterization of the sample was giyen in 
IZarattini et al.l (I20l4 hereafter FOGO IV). The structure of this 
paper is as follows. Section |2]is deyoted to describing the sam¬ 
ple; Sect.|3]shows how the LFs are calculated; Sect. |4]describes 
the dependence of the LFs on the magnitude gap; and Sects. |5] 
andl^present the discussion and the conclusions, respectiyely. 

For this work, the adopted cosmology is Hq - 70 km s ' 
Mpc *, Oa = 0.7, and Q.m - 0.3. 
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2. Description of the sampie 


We used two samples of galaxy aggregations to analyze their 
LF and their dependence on the magnitude gap. The first sam¬ 
ple (hereafter SI) is composed of 34 g roups and clusters of 
galaxies selected by ISantos et al.l (1200711 and analyzed in de¬ 
tail in FOGO IV. These systems were selected as FG candidates 
from the Sloan Digital Sky Surv ey Data Release 5 (SDSS DR5; 
lAdelman-McCarthv et al.ll2007h . and they present a wide range 
in redshift (0 < z < 0.5). However, in our detailed analysis of 
this sample using deep r-band images and multi-object spec¬ 
troscopy, it is shown that only 15/34 system s meet at least one 
of t he definitions of fos sil systems given by iJones et ^ (l2003h 
and iDariush et al.l(l2010l) . The former authors claim that a cluster 
or group of galaxies is fossil if it has a gap in magnitude larger 
than 2 in the r-band between the two brightest member galaxies 
(Ami 2 > 2). The latter authors suggest that a system is fossil if 
the gap in the r-band between the first- and fourth-ranked galax¬ 
ies is larger than 2.5 (Am^ > 2.5). Both quantities are defined 
within half the virial radius. 


The need for a second sample comes from the mean value of 
Ami 2 being ~ 1.5 in the SI sample. Only four systems have gaps 
lower than 0.5. We are interested in studying the dependence of 
the LF on the Amn key parameter, so we need to extend the 
sample toward systems with small Aoti 2 . 

For these reasons, we used a second sample (hereafter 
S2), taken by lAguerri et alJ (l2007h . These systems were se¬ 
lected as all the galaxy agg regations wi t h known redshift at 
z < (3.1 from the catalo g s of Abell et ^ (Il989h, Zwdc W et alJ 
(1196Ih . iBohringer et alJ (l2000h . and \^ges_etaL ( 1999|) that 
were mapped in the SDSS DR4 (lAdelman-McCarthv et alJ 
l2006l) . This selection results in 88 systems. Of these, we used 
only those for which the two brightest members were spectro¬ 
scopically confirmed, which limits the final number of systems 
from this sample to 76. For this sample, the mean value of Ami 2 
is ~ 0.7, with a standard deviation of ~ 0.5. 


The general prop erties of the two samples ar e presented in 
IZarattini et al.l (12014l) for the S1 sample and in lAguerri et al.l 
( 2007ll for the S2 sample. Both samples are map ped in the SPSS, 
and w e used their r-band model magnitude (see lStoughton et ^ 
I 2 OO 2 I) . The selection was done using SDSS DR5 for the SI 
sample and DR4 for the S2 sample, but the magnitude used 
for this work were ta ken from the more recent SDSS-DR7 
(lAbazaiian et al.l l2009t) . These magnitudes were corrected for 
galactic extinction and K-correction. The former was obtained 
by using the r-band extinction par ameter provided by SDSS. The 
latter was computed follow i ng thelChilingarian et al.l (1201 Ol) and 
IChilingarian & ZolotukhinI (12012l) prescriptions. Moreover, the 
R 200 radius of each system was computed using X-ray data from 
the ROSAT satellite (see FOGO Ill and FOGO IV). The LFs of 
this work were computed within half the obtained R 200 radius. 

The X-ray luminosity of SI ranges between 41.9 < 

log (Lx.si [erg s“']) < 45.1, whereas that of the S2 sample 
varies in the range 41.6 < log(Lx,S 2 [ergs^*]) < 45.2. The me¬ 
dian values of the X-ray luminosity of the SI and S2 samples 
are log(Lx,si [ergs^*]) = 44.1 + 0.7 and log(Lx..s 7 [erg s~‘)l ^ 
43.8 + 0.7. Masses can be obtained using Eq. 3 of iRvkoff et al.l 
(l2008l) . after an adequate cosmology correction, and it varies 
in the ranges 13.1 < log(M 2 oo,si [Mq]) < 15.0 and 12.9 < 
log (M200,S2 [Mq]) < 15.1 for the SI and S2 samples, respec¬ 
tively. The median values of the masses are log (M 2 oo,si [M©]) = 
14.6 + 0.4 and log (M 2 oo,S 2 [Mq]) = 14.4 + 0.4, respectively. 
Finally, the velocity dispersion of the SI and S2 galaxy sys¬ 
tems span the ranges 250 < cr,, si [kms-‘] < 1200 and 250 < 



(mag) 


Fig. 1. Galaxy background estimations. Red filled circles rep¬ 
resent the background used in this work, black asterisks are 
taken f rom lCanak et al.l (120041), blue triangles f rom fYasuda et al.l 
(l200lh . vi olet stars fromi Huung et"^ (1200ih . and green plus 
signs from iMetcalfe et akl (1200 ll) . 


o‘v,s 2 [kms '] < 1000, respectively. The median values are 
759 + 253 km s“* and 557 + 170 km s“', respectively. The 
Kolmogorov-Smirnovtest was applied to the X-ray luminosities, 
masses, and velocity dispersion. Results from this test indicate 
that the SI and S2 samples do not come from the same parent 
distribution. However, we are not directly comparing these two 
subsamples, since we want to compute the LFs in bins of Anin- 
In Sect. 15.II we show that the four analyzed subsamples actually 
come from the same parent distribution. 

2.1. Magnitude gap determination 

For determining Am^ in the S1 sample we proceeded as follows. 
We considered the four brightest galaxies of the systems within 
the spectroscopically confirmed members and possible members 
(see FOGO IV for both definitions), and then for each of these 
galaxies, we computed the magnitude as the mean value of three 
different magnitudes. W e used the model and th e Petrosian mag¬ 
nitude from the SDSS (IStoughton et akl l2002l) and the MAG- 
BEST magnitude by ob tained analyzing our ow n deep r-band 
images with SExtractor (iBertin & Arnouts|ll996f) . We used this 
mean value to compute both Aoti 2 and A»ji 4 and used the stan¬ 
dard deviation of the mean value for computing the uncertainties. 
The detailed procedure can be found in the EOGO IV paper. 

Eor the S2 sample, we used the same methodology to ob¬ 
tain Ami 2 and Ami 4 except that we did not have our own photo¬ 
metric images. Thus, we used the mean value of the model and 
Petrosian magnitudes alone. 


3. Galaxy luminosity function determination 

3.1. Luminosity functions of individuai systems 

There are two methods that are widely used in the literature to 
compute the LEs of individual systems: the spectroscopic and 
the photometric ones. The former is, in principle, the most accu¬ 
rate. It is based on an extended redshift catalog, which allows for 
a detailed study of the system membership. Nevertheless, it re¬ 
quires a large amount of observational time and generally larger 
telescope apertures than the photometric one. This method has 
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mainl y been applied to the study of nearby i ndividual clusters 
(e.g., [Rines & Gellej20()^lAgulli et al.l2()14ll . The photometric 
method requires less observing time and consists in computing 
the galaxy number counts, as a function of magnitude, within the 
system region and in a region of the sky in which no structure 
is present. The difference between these two quantities repre¬ 
sents the galaxy system LF. It is a statistical method, so its main 
problem is that wherever the background field is located, it will 
not be exactly the same background of the system itself, mainly 
because of cosmic variance. Moreover, if the system is poor, it 
could have very limited contrast with respect to the background, 
leading to large uncertainties in the LF. 

The LFs of individual systems, computed in half the R 200 
radius, were obtained by using a hybrid method. This procedure 
uses both the photometric and spectroscopic information that we 
have for each system. The galaxy LF of each system in the j-th 
magnitude bin is given by 

cPj^N^j + iNj-N,j)xPj, ( 1 ) 

where N^j is the number of spectroscopically confirmed mem¬ 
bers, Nj the total number of galaxies, and N^j the number of 
galaxies with recession velocity measurements. Finally, Pj is 
given by 

Pj = (Nj - Nbj)/Nj, (2) 

where Nb.j is the number of background galaxies. We notice that 
Pj represents the probability that a galaxy would be considered 
in the LF when only photometry is available. The background 
was obt ained by averaging four different fields u s ed in the lit¬ 
erature d Yasuda et al.| |200lt iMetcalfe et al.l 1200 it iHuang et al.l 
l2001t ICanak et alJl2004l) . The resulting background distribution 
is shown in Fig. [1] The completeness limit of our LFs is r = 21.5 
mag, and it is a conservative choice since the nominal complete¬ 
ness of SDSS DR7 r-band is 22.2 mag. 

The advantage of using this methodology is that we obtained 
a quasi-spectroscopic LF for the brightest bins, where the mag¬ 
nitude gap arises. In fact, the SI sample is ~85% complete down 
to m^ = 17 and the S2 sample i s ~9Q% complete do wn to mr=17.5 
(see F ig. 4 of IZarattini et alJ (1201 4l) and Fig. 1 of lAguerri et alJ 
( 2007h '). Then, when we move to fainter magnitudes, the number 
of measured redshifts decays rapidly and the LF is dominated 
by the statistical background subtraction. What we obtained is a 
hybrid LF that is more accurate at the bright end than the photo¬ 
metric one, but that is not as time consuming in terms of observa¬ 
tions as a full spectroscopic LF in the faint end. The uncertainty 
associated to the LF is calculated using the error propagation of 
the terms in Eq. [Hand adding in quadrature the cosmic variance 
following lHuan^t al.l ([1993- 

In Fig.l^we show some examples of the individual LFs cal¬ 
culated with this procedure. For some rich objects, the individ¬ 
ual LFs are clearly defined (e.g., Abell 1066 or FGS02), but for 
poorer systems the LFs present large uncertainties (e.g., Abell 
724 or FGS15). All the LFs presented in this work were com¬ 
puted excluding the BGGs, as is usually proposed in the litera¬ 
ture. 

3.2. Stacked luminosity functions 

We stacked the individual LFs for all the SI and S2 systems 
with z < 0.25 to deal with small numbers. This cut in redshift 
was needed to guarantee that all systems will reach at least a 
magnitude of Mr - -19.5 (given a completeness limit of the 
SDSS of nir - 21.5) so that they can be normalized using the 


region -21 < Mr < -19.5. As noted by other authors (e.g., 
iPonesso et al.ll2005h . stacked LFs are not only a useful tool for 
checking the universality of the LF, but they are useful for calcu¬ 
lating the LF of systems with high accuracy when the individual 
ones have poor statistics. 

There are different methods i n the lite r ature o f stacking LFs. 
We used the one proposed by ICollessI (Il989l) . in which the 
stacked LF can be obtained by combining the individual ones 
according to the formula 


_ ■y ^ 
“ m, 4^ Wo’ 


(3) 


where is the number of galaxies in the j-th bin of the stacked 
LF, Nij the number of galaxies in the j-th bin of the i-th individual 
system’s LF, Wo the normalization of the i-th system LF in the 
region -21 < Mr < -19.5, my the number of systems contribut¬ 
ing to the j-th bin, and Wo is the sum of all the normalizations 

(Wo = li/Wo)- 

The formal errors are computed according to 



where S(pcj and 6Nij are the formal errors in the LF’s j-th bin for 
the composite and i-th system, respectively. 

The Colless method can be safely used under a few condi¬ 
tions; first of all, the magnitude limit of all stacked systems must 
be fainter than the values used for normalization (Mr - -19.5 
in our case). Second, the normalization region should be large 
enough to be representative of the richness of the systems. 
Finally, the number of systems contributing to each bin should 
be statistically relevant. 

As introduced at the beginning of this section, we applied 
this method to all the systems in our sample with z < 0.25. Using 
this subsample, we satisfied all the requirements for applying 
the Colless method. The total number of systems turns out to 
be 102, and the resulting composite LF is presented in Fig. [3] 
We fit neither a single nor a double Schechter function to the 
data, because the Spearman test told us that none of them was 
representative of the data. We fit an exponential function to the 
faint end of the LFs using the last five points. The form of the 
exponential is 10^'", where m represents the magnitude and k is 
related to Schechter’s a parameter by 



(5) 


The resulting faint-end slope is a = -1.27 + 0.11. Hereafter, 
all the presented exponential slopes were fitted using the five 
faintest points of each LF. 


4. Dependence of the luminosity function on the 
magnitude gap 

We divided the sample of the 102 systems with z < 0.25 into four 
subsamples, which differ from one another in the value of Ami 2 . 
The first subsample is composed of 31 systems with Ami 2 < 0.5, 
the second of 24 systems with 0.5 < A»ji 2 < 1, the third of 26 
systems with 1 < Amn < 1.5, and the fourth of 21 systems with 
Aoti 2 > 1.5. This division is arbitrary and was done in order to 
have a statistically significant number of systems in each range 
of Ami 2 and to trace the dependence of the LF with Ami 2 in the 
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-26 -24 -22 -20 -18 -16 -26 -24 -22 -20 -18 -16 -26 -24 -22 -20 -18 -16 -14 

Mr (mag) 


Fig. 2. Examples of individual LFs for 9 systems taken from our sample. For the more massive ones, such as ABEFF0085 and 
FGS02, the FFs are determined well, whereas for less massive systems, such as ABEFF0724 and FGS15, the FFs have large 
uncertainties. 


Table 1. Best-fitting parameters of a Schechter fit to the regular 
(top) and relative (bottom) FFs. Reported uncertainties represent 
the 99% confidence level (c.l.) of each parameter. 
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best possible way. In Fig.|4]we show the stacked FFs for the four 
subsamples. Qualitatively, the slope of the bright end is similar 
for the four FFs, but not the faint-end one. In particular, the sys¬ 


tems with Ami 2 < 0.5 show a steeper faint end than those with 
Ami 2 > 1.5, whereas the other two subsamples represent inter¬ 
mediate cases. To quantify this effect, we fit a single Schechter 
function to each FF shown in Fig.H] In this case, the Spearman 
test confirms that a single Schechter function is a reasonable rep¬ 
resentation of the data of eac h subsampl e. _ 

The Schechter function (ISchechteijll976ll is the most ac¬ 
cepted expression to describe the galaxy FF parametrically. Its 
formulation can be written as follows: 

where (ff is a normalization factor defining the overall density 
of galaxies and M* is the characteristic magnitude. The param¬ 
eter a describes the faint-end slope of the FF, and it is typically 
negative. The results of the fit are shown in the upper part of 
Table [1] There are differences in both the bright and the faint 
ends. In the former, the larger the gap, the fainter the M*. In the 
latter, the larger the gap, the flatter the a. We plotted the 68%, 
95%, and 99% confidence level (c.l.) contours for M* and a of 
the Schechter fit in Fig.|5] We refer to FFs computed using this 
method as “regular” FFs. 

Moreover, for the faint end, we fit an exponential function 
as well. This check is useful for two reasons: it helps to quantify 
the effect of the known degeneracy between M* and a, and it can 
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be useful to compare the results with other LFs that are not well 
described by a Schechter function. The obtained values for the 
exponential faint-end slope are a = -1.43+ 0.12, -1.39 + 0.17, 
-1.18+0.31, and -1.02+0.32 for the four LFs with an increasing 
gap, respectively. These values are higher than for Schechter’s a 
parameter, but the trend is the same. However, the differences in 
a between the four slopes are not statistically significant because 
of the large uncertainties of the stacked galaxy LFs for systems 
with larger Ami 2 . 

The BGGs in ou r sample span a four-magnitude range (see 
IZarattini et alj|20f^ . This can affect the shape of the bright end 
of the stacked LFs, because the individual LFs are not aligned in 
magnitude. To avoid any effect associated with the stacking of 
different galaxy populations, we have computed the stacked LFs 
as a function of the relative magnitude. This was obtained by 
calculating the differences between the magnitude of the galax¬ 
ies and the BGG (AMr - Mr - Mr^Bcc) for each system. Using 
this method, all the BGGs are located at AMr - 0. We refer to 
the resulting LF as the “relative” LF. In this picture the value of 
M* loses its physical meaning, but we are interested in highlight¬ 
ing the differences in the Schechter parameters between the four 
subsamples, not in their absolute values. 

There are two differences between the methodology that we 
applied for the regular LF and for the relative LF. The first dif¬ 
ference is the already mentioned shift of the magnitudes in the 
calculation of relative LFs. The second difference is that, for the 
relative LF, spectroscopically confirmed zero-galaxy bins enter 
into the stacking formula (Eq. [3]) for each bin of magnitude and 
for each system. In contrast, for the regular stacking, this was not 
possible because of the large difference in the magnitude of the 
central galaxies. In fact, the stacking is actually a computation of 
the mean value of the LF in each bin, but when the dominant pa¬ 
rameter is the absolute magnitude, mixing massive systems with 
small groups could be misleading. If a group is dominated by 
a central galaxy of = -22 and has, for example, a spectro¬ 
scopically confirmed two-magnitude gap, when we try to stack 
it with a massive cluster whose central galaxy has = -25, the 
spectroscopically confirmed gap of the group would affect a part 
of the cluster that is three to five magnitudes fainter than the cen¬ 
tral galaxy of the cluster itself. This part, assuming a Schechter 
profile for the distribution of galaxies, would probably be lo¬ 
cated beyond the elbow (M*) of the LF of the cluster. Thus, we 
expect that in this region the cluster presents a large number of 
objects. Doing the stacking in this case would imply reducing 
the galaxies that are present in that bin by a factor of 2. Clearly, 
if the number of systems is more than two, as in our case, the 
effect would be softened, but we expect a flattening of the elbow 
region if we do not take this aspect into account. 

In Fig.|6]we show the stacked relative LF for the whole sam¬ 
ple of 102 systems with z < 0.25. As we did for the regular LF, 
we fit an exponential to the faint end of the relative one. The re¬ 
sulting slope isa = -1.25+0.09, which is compatible (within the 
uncertainties) with the value measured for the regular LF. Figure 
|7] shows the relative LFs of the four subsamples with different 
magnitude gaps. Qualitatively, we found differences in both the 
bright and the faint ends of the four stacked relative LFs. We fit 
a single Schechter function to these relative LFs and show the 
obtained M* and a parameters in Table [T| Their uncertainties 
are shown in Fig. [8] The Schechter parameters of the smallest 
and largest magnitude gap regimes have a greater difference in 
value for the relative LFs than the regular LFs. Once again, we 
fit an exponential to the faint end of the four LFs. The resulting 
faint-end slopes are -1.37 + 0.12, -1.43 + 0.13, -1.24 + 0.22, 
and -0.95 + 0.17, moving from the smaller to the larger gap. In 


this case, the exponential fits also show differences at the 1-cr 
level between the systems with < 0.5 and Amn > 1.5. The 
larger uncertainties are due to the small number of data points 
available for the exponential fit. This fit has been performed as a 
test to break the degeneracy between M* and a that is obtained 
when both the bright and faint ends are fitted at the same time. 
The exponential fit of the LF faint end is statistically less signif¬ 
icant. But, it has to be considered as a double check of the slope 
derived by fitting the full LF with a Schechter function. 

In Fig. 121 there are some points in the LFs that are located at a 
magnitude difference from the BGG that is smaller than the mag¬ 
nitude gap used to dehne the four subsamples. For example, in 
the sample with Aoti 2 >1.5 the LF starts at Aoti 2 =L This appar¬ 
ent contradiction is due to the definition of the magnitude gap. 
It is obtained by only using possible members. The selection of 
possible members cannot be made using a statistical background 
subtraction, because we need individual galaxy information. For 
this reason, in the determination of the gap we used a gener¬ 
ous cut in ph otometric redshift to account for possible cluster 
members (see IZarattini et al.]|20l4 for details). Nevertheless, no 
photometric-redshift cut was applied for the computation of the 
galaxy LF. The methods are not incoherent, since the error bars 
of these peculiar points are always compatible with zero. These 
points should be considered as statistical fluctuations in the very 
bright part, but they do not affect the results owing to their large 
uncertainties. 

4.1. Dwarf-to-giant galaxy ratio 

It is important to remember that the stacked LFs is arbitrarily 
normalized, which means that looking at the absolute number of 
galaxies at both the bright and the faint ends of the LFs can be 
misleading. For this reason, we analyzed the so-called dwarf-to- 
giant galaxy ratio (DGR) for the four subsamples. We defined 
as giant galaxies those in the range -22.5 < Mr < -20, and 
as dwarf galaxies those in the range -19 < Mr < -17. The re¬ 
sulting DGR for the four subsamples with increasing magnitude 
gaps are 2.27 + 0.03, 1.70 + 0.02, 1.37 + 0.02, and 1.33 + 0.02. 
The DG R is very diffic u lt to c ompare with the literature. For 
example, iPonesso et alJ (l2005l) defined galaxies in the range 
-18 ^ Mr < -16.5 as dwarfs and those with Mr < -20 as gi¬ 
ants, ISanchez-Janssen et alJ (l2008l) considered as dwarfs those 
with Mr > Ml + 1 an d as giants those with Mr < M* , and 
IWeinmann et al.l (1201 ill defined as dwarf galaxies those with 
-16.7 > Mr > -19 and as giant those galaxies with Mr < -19. 
While the exact value of the DGR is therefore highly arbitrary, 
the important result of this exercise is to show that we recover 
the trend previously found using the fits to the LF; the relative 
number of dwarfs systematically decreases in systems with pro¬ 
gressively larger magnitude gaps. 

5. Discussion 

5.1. Caveats of the results 

In Figs.|4]and|7]we have stacked all the available systems, mixing 
clusters, and groups. This can affect the result, since in some of 
the four magnitude-gap bins we could be dominated by massive 
clusters, while in others the dominant systems could be groups. 
Differences can be found in the literatur e between the Schechter 
param eters of clusters and groups (e.g., IZandivarez & MartmezI 
1201 ll) . To test this aspect, we ran a Kolmogorov-Smirnov test, 
which confirmed that the distributions of <t^ (which is a mass 
proxy, as suggested bv iMunari et al.l (l2013h . and it has been ob- 
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Fig. 3. Stacked LF for all the systems of our sample with z < 
0.25. The solid lines represent the faint-end slope for a - -2.0, 
which is the value predicted by standard CDM theories, a - 
-1.3, obtained from our fit, and a = -1.0, which is the value for 
a flat LF. 

tained from Lx) for the four subsamples that are not different. 
Moreover, we computed the median o-y for the four subsam¬ 
ples. The resulting values are o-y = 557, 591, 587, and 545 km 
s ', with standard deviations of 171, 159, 206, and 200 km s ', 
for the four subsamples ordered with increasing magnitude gap. 
These two tests indicate that the four subsamples show the same 
mean velocity dispersion, hence the same mean mass. We also 
divided the whole sample (S1+S2) into two mass bins, defining 
as groups those systems with o-y < 580 km s“' and as clusters 
those with o-y > 580 km s '. This value represents the median 
value of o-y calculated over the whole S1 -i-S2 sample. In Fig. |9] 
we present the stacked LFs according to the mass of the sys¬ 
tems and the ratio between the two LFs. It can be seen that the 
LFs are in good agreement with one another, because the dif¬ 
ferences are always compatible with zero except for one point, 
M,. = -18.5. In fact, in this point there is a dip in the more mas- 
sive systems, as already fou nd by other authors (see Fig. 11 in 
iTrentham & Hodgkinll2002h . In conclusion, the observed differ¬ 
ences in the LFs do not seem to be related to the mass distribu¬ 
tion of the systems in the four subsamples. 

As pointed out in Sect. 13.11 we used a galaxy background in 
order to compute the galaxy LFs. We tested how the results were 
affected by changing the adopted galaxy background. We used 
the 26 systems of the SI sample that are closer than z = 0.25. 
These systems are representative of the whole S1+S2 sample in 
terms of mass (o-y), and their /?200 are smaller than 15 arcmin. 
We then calculated a local background for each system in an an¬ 
nulus between two to four times their 7?200 radius. We divided 
the annulus into 20 regions of the same area, as proposed by 
iPonesso et al.l (l200-5h . Then, we counted the number of galaxies 
for each bin of magnitude in each sector, and Anally we calcu¬ 
lated the mean value of the local background by averaging all the 
sectors, using a sigma-clipping algorithm to exclude sectors that 
are at more than 3-cr from the global mean value. In Fig.fTOlwe 
show the values of the local background for each system and the 
1-cr and 3-0- uncertainties of the global background presented in 
Fig. [T] It can be seen that up to m,. = 17 the local and global 
backgrounds seem to disagree with one another, but this does 



Mr (mag) 


Fig. 4. Lower panel; Stacked LFs of systems with different gaps 
in magnitude. Black open squares are systems with gap < 0.5, 
red filled squares represent systems with 0.5 < Am 12 < 1.0, 
violet open circles indicate systems with 1 < Ami 2 < 1.5, and 
green filled circles are systems with gap >1.5. The four LFs have 
been moved by an arbitrary offset for display purposes. Upper 
panel; histogram of the number of systems that are contributing 
to each bin. The color code is the same as in the lower panel. 

not present a problem in our case because, as we already men¬ 
tioned, we are using a hybrid method for the LF, and both the S1 
and the S2 samples have a spectroscopic completeness of more 
than 85% up to m^ = 17. Moreover, for magnitudes brighter 
than mr -17, the local background method is less reliable, since 
bright galaxies are scarce, and a large area is needed to properly 
take them into account. For magnitudes between -17 and 
ftir = 21.5, which is our conservative completeness limit, the 
global and local backgrounds are in good agreement, so no large 
differences are expected in the calculation of the LF by varying 
the background. 

The stacking procedure computes a mean LF, using for each 
bin of magnitude only those clusters where the magnitude limit 
is fainter than each specific bin. Thus, each bin of the compos¬ 
ite LF is formed by a different number of averaged points. At 
the faint end, the systems that are contributing are the closest in 
terms of redshift, since our sample is limited in apparent mag¬ 
nitude. However, the faint-end slop e is not expected to c hange 
in the redshift range 0 < z < 0.25 (iGozaliasl et al.ll^l4tl . thus 
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Fig. 5. Uncertainty contours for the Schechter fits of Fig. |4] 
Contours represent 68%, 95%, and 99% c.l. and the color and 
symbol codes are the same as in Fig. |4] The error bars are the 
1-cr uncertainties of the Schechter fit as reported in Table [1] 



Fig. 6. Stacked relative LF for all the systems of our sample with 
z < 0.25. The two solid lines represent the same a values as in 
Fig. [3] 

we think that the method does not introduce any bias into the 
dependence of the faint-end slope on the magnitude gap. In the 
bright end, differences could arise due to the use of our hybrid 
method. Nevertheless, no differences are expected in the bright 
end when applying a fully photometric or a fully spectroscopic 
method, although the latter has smaller uncertainties. Thus, the 
use of a hybrid method should not affect the computation of the 
M* values, but it should help in reducing their uncertainties. We 
have plotted in Figs.|4]and|7]the histograms showing the number 
of systems per bin of magnitude used in computing the stacked 
LFs. 

We also analyzed the differences in the redshift distribu¬ 
tions of the four subsamples. The Kolmogorov-Smirnov test in¬ 
dicates that the four redshift distributions differ from one an¬ 
other. The median redshift value for subsamples with increasing 
magnitude gap are z - 0.064, 0.077, 0.088, and 0.11, respec¬ 
tively. Nevertheless, the differences in redshift are small, and 



Fig. 7. Same as Fig.|4]but with relative LFs. 


iGozaliasl et al.l (12014li show that no evolution is expected in the 
faint-end slope for both fossil and non-fossil systems since z = 1. 
Moreover, the lookback time at z = 0.1 is ~ 1 Gyr, which is a 
small amount of time to see an evolution. Thus, we think that the 
differences in redshift in the four subsamples do not represent a 
bias for our results. 

To conclude our analysis of the possible caveats, we investi¬ 
gated how the uncertainties in the magnitude determination de¬ 
pend on the magnitude itself, and how this can affect the com¬ 
putation of the LFs. We used SDSS model magnitudes, so to 
constrain the uncertainties we analyzed the distribution of the 
modelMagEtTr parameter. The median uncertainty at m^ = 21.5 
- our conservative completeness limit - is 0.15 mag. Since our 
bins are 0.5 mag wide, we expect that the photometric uncertain¬ 
ties do not affect the results. 

5.2. Comparison with the iiterature 

We can use the regular LFs to compare the results of this work 
wit h other results in the literature for the r-band. For exam¬ 
ple, [Popeiio^eral] (l2005h found a slope of a = -1.30 + 0.06 
for the bright part (M^ < -18) of the stacked LF of ~100 
clusters inside 1 Mpc, and a slope of a = -1.29 + 0.09 for 
the sa me sample within 0.5 Mpc. In contrast, Ide Filippis et al.l 
(1201 ih find a faint-end slope of a = -0.99gg2 for ^ stacked 
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(mag) 

Fig. 8. Uncertainty contours for the Schechter fits of Fig. |7] 
Contours represent 68%, 95%, and 99% c.l. and the color and 
symbol codes are the same as in Fig. |4] The error bars are the 
1-cr uncertainties of the Schechter fit as reported in Table [1] 



Mr (mag) 


Fig. 9. Stacked LFs for both clusters (black filled rectangles) and 
groups (red open rectangles) in the upper panel. Here we consid¬ 
ered objects with cTy < 580 km s“' as groups and objects with 
(Tv > 580 km s ' as clusters.The ratio of the two LFs is plotted 
in the lower one. 



Fig. 10. Local background (black dots) superimposed on the l-cr 
and 3-cr contours of the global background (red and gray shaded 
areas, respectively). The vertical dotted line is the completeness 
limit of this work. 

LF of ~1500 systems. We can also try to compare our stacked 
LF with individual systeni s analyzed in the literature. For exam¬ 
ple, [WneF&Geii^ yOOSl) analyzed spectroscopic LFs of Abell 
2199 and the Virgo cluster, finding a — -1.02 + 0.05 for the 
former and a - -1.28 + 0.06 for the latter. These results, ob¬ 
tained with different techniques, are compatible with our faint- 
end slo£e_of_Q_=_-L27± 0.11, wi th the only exceptio n be¬ 
ing Ide Filippis et al.l (1201 Ih . Finally, iBarrena et al.l (120121) find 
no difference in the a parameters in relaxed and unrelaxed clus¬ 
ters with a - -0.86 + 0.27 and a = -0.99 + 0.21. This result 
is of particular interest since fossil systems are thought to be old 
and, thus, more relaxed than non-fossils. 

The evolutionary state of the system could be seen as an 
explan ation of our results, too. In fact, Ilglesias-Paramo et al.l 
(l2003h studied two clusters (Coma and Abell 1367) with sim¬ 
ilar X-ray luminosities and redshifts and found significant dif¬ 
ferences in their faint-end slopes. The authors suggest that these 
differences could be explained by differences in the evolutionary 
state of the two clusters. Thus, if applying this result to our case, 
a possible explanation for the observed differences is that the 
Ami 2 parameter could be an indicator of the evolutionary state 
of a system. Either way, extended work on the substructure in 
fossil systems remains to be done, so differences in their evolu¬ 
tionary state cannot be proved. 

It is also interesting to compare the result of the Ami 2 >1.5 
subsample with the results found for fossil systems in the liter¬ 
ature. In fact, in this subsample, 13 out of 21 systems are spec- 
troscopically confirmed fossils. iKhosroshahi et al.l (l2006h find 
a = -0.61 + 0.20 for t he background-corrected photom etric 
LF within 0.5 R 200 and iMendes de Oliveira et al.l (l2006l) find 
a - -0.64 + 0.30 in ~0.3 R 2 Q 0 , whereas in FOGO 1 we found 
a = -0.54 + 0.18 for LFs within 1 Mpc. It can be seen that, 
despite the large uncertainties, all the results seem to point to a 
value of the faint-end slope that is higher than -1 for fossils. In 
this sense, our result agrees with the literature, since our faint- 
end slope in 0.5 R 20 Q for systems with large magnitude gaps is 
a - -0.78 + 0.12 for the regular LF and a - -0.77 + 0.14 
for the relative LF. That our result is the lowest one can be inter¬ 
preted as an effect of non-fossil systems being in this subsample, 
owing to its definition. The effect of this contamination would 
be a steepening of the LF. This agrees with our general result 
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that the faint-e nd slope increases with the gap in magnitude. In 
a recent work. lLieder et al.l (1201 3h analyzed the FG NGC 6482. 
They present a deep spectroscopic LF, down to Mr — -10.5. The 
faint-end slope they found was a - -1.32 + 0.05 in one virial 
radius, which is steeper than the other works in the literature. 
However, their result is not directly comparable to the others, 
since they fit the faint-end slope do wn to a much deepe r magni¬ 
tude. Finally, in a very recent work. lWen & HanI (1201 5h have an¬ 
alyzed the dependence of the bright end of the LF o n the cluster 
dynam ical state. They used the method presented in I Wen & HanI 
(1201 3h to create three subsamples of clusters with different dy¬ 
namical states. They conclude that more relaxed clusters have 
fainter M*. Thus, our results for the bright end would indicate 
that systems with a larger magnitude gap would be more dy¬ 
namically relaxed. Nevertheless, a general study of the dynami¬ 
cal state as a function of the magnitude gap remains to be done. 

It is worth noticing that the LFs of the four subsamples— 
divided for different magnitude-gap bins and both in the reg¬ 
ular and relative cases—are reasonably represented by a sin¬ 
gle Schechter function. Nevertheless, the stacked regular and 
relative LFs of the 102 systems are not well fitted by a sin¬ 
gle Schechter function. These differences cannot be due to the 
method, since only stacked LFs are analyzed. Thus, we interpret 
this result as another hint that the observed differences in the M* 
and a parameters of the four subsamples are real and not due to 
statistic effects. 

5 . 3 . Implications for formation scenarios 

The main results of this work suggest that there is clear depen¬ 
dence on the magnitude gap in the bright part of the LF and a 
less significant dependence on the magnitude gap in the faint- 
end slope of the LF. The characteristic magnitude of the LF can 
be interpreted as the mean luminosity of the bright galaxy pop- 
ulation of a system, when the BGG is excluded from the LF 
(ICooray & Milosaylieyidl2005h . Thus, the obseryed difference 
of 1.3 magnitudes in the M* parameter of the regular LF cor¬ 
responds to a factor ~3 in the mean luminosity. These results 
are consistent with the most accepted scenario for the formation 
of the magnitude gap, namely that it deyeloped through dynam¬ 
ical friction, and a ll the M* galaxies merge d in a single, mas- 
siye central object (ID’Onghia & Lakell2004li . Howeyer, it is not 
clear whether the magnitude gap of these systems was form ed 
at hi gh redshift o r more recently (iDfaz-Gimenez et ^l2008h . In 
fact. iRaouf et al.l (1201 4h suggest that Am 12 alone is not a good 
age indicator. They claim that there is a trend in age from groups 
to clusters, where the former are, on ayerage, older than the lat¬ 
ter. 

Using their dating method, our sample of fossi ls would be 
mainly dominated by young systems. Moreoyer, ISmith et^ 
(1201 oh suggest that the formation of the large magnitude gap 
could depend on both the formation time or the recent infall his¬ 
tory of the systems. Thus, a fossil system could also form in a 
recent epoch, and it could eyolye into a regular system in the 
future, by inte racting with other groups/clusters in the same re- 
gion (see also |yon Benda-Beckmann et al.ll20'0^ . Neyertheless, 
iDeason et al.l (1201 3h show that, on ayerage, older and more con¬ 
centrated halos haye larger mass (magnitude) gaps, but the scat¬ 
ter is important because of the transient nature of the satellite 
population. 

These two results show that the magnitude gap cannot be 
used as an unambiguous age estimator for indiyidual systems, 
but that there is a statistical trend in the magnitude gap-age re¬ 
lation. In this work we circumyent this limitation by comput¬ 


ing stacked LFs for systems with similar magnitude gaps. For 
this reason, it is reasonable to assume that the systems in the 
larger Ami 2 bin are, on ayerag e, dynamically older tha n the oth¬ 
ers. But not only age matters. ISommer-LarsenI (l2006h suggests 
(using smoothed particle hydrodynamics simulations) that the 
anisotropy of the orbits can play a role as well. This kind of orbit 
can bring the infalling population close to the center of the po¬ 
tential well, thus fayoring the merging of massiye galaxies. They 
conclude that the more radial the orbits at the time of formation, 
the more fossil the system will be at the present time. 

The obseryed dependence of the faint-end slope of the LF 
on the magnitude gap is more intriguing. There is some eyi- 
dence in the literature showing that some nearby galaxy clusters 
(e.g., the Coma and the Virgo clusters) conta in a small number of 
dwarf galaxies in th e innermost regions ('see lAguerri et al.ll2004l : 
iTruiillo et alj|2002h . Neyertheless, our results point out that the 
number of dwarfs in the innermost regions (R < 0.5 R 200 ) de¬ 
pends on the magnitude gap as well. The low masses of dwarf 
galaxies make them less susceptible to dynamical friction. In ad¬ 
dition, the large yelocity dispersion in galaxy clusters makes the 
merging of dwarfs in the recent epoch a yery rare eyent. 

One appealing possibility is that the paucity of dwarfs in sys- 
te ms with large gap s is rel ated to the more radial orbits predicted 
by lSommer-LarsenI (l2006h . These eccentric orbits can efficiently 
bring infalling groups close to the center of the potential well. 
Once there, not only will the more massiye galaxies of these 
groups merge on a relatiyely short timescale, but strong tidal 
forces can efficiently disrupt lower mass halos. Furthermore, 
if the assembly occurs at relatiyely early times, the suryiying 
subhalos will spend a significant amount of their history orbit¬ 
ing within rather massiye halos—thus increasing the chances 
of eyentual disruption. A si milar tidal force-driyen disruption 
mechanism was proposed by lLonez-Cruz et ^ (Il997h to explain 
the flattening of the faint-end slope of the LF in the central re¬ 
gions of clusters hosting cD centrals. 


6. Conclusions 

We analyzed a sample of 102 systems with redshift z < 0.25 
in order to determine the properties of their LFs. The sample 
was diyided into four subsamples, coyering different ranges of 
Ami 2 . In particular, the first subsample included systems with 
Ami 2 < 0.5, the second systems with 0.5 < Ami 2 < 1.0, the 
third systems with 1 < Ami 2 < 1.5, and the fourth systems 
with Am 12 > 1.5. The LFs were computed in half the R 200 ra¬ 
dius using a hybrid method, which allowed us to use both pho¬ 
tometric and spectroscopic data. Moreoyer, to better define the 
differences in the parameters between the four subsamples, we 
calculated the relatiye LF. For each galaxy system, the relatiye 
LF was obtained by shifting the LF in magnitude such that the 
BGG magnitude is zero. 

Our results can be summarized as follows: 

- The faint-end slope of the regular stacked LFs of the 102 
systems in our sample turns out to be cr = -1.27 + 0.11. The 
slope of the relatiye stacked LF is or = -1.25 + 0.09. The 
two slopes are in good agreement. These slopes haye been 
obtained by fitting an exponential function and are calculated 
using the last flye points of each LF (-18 < Mr < -16.5 for 
the regular LF and 4.5 < AMr < 1 for the relatiye LF). We 
fit an exponential to these stacked LFs because they are not 
adequately described by either a single or a double Schechter 
function. 
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- In both the regular and relative stacked LFs, if we divide 
the sample of 102 systems into four subsamples of different 
Ami 2 , the M* values of the Schechter fit change. In partic¬ 
ular, the larger the gap, the fainter the M*, as is expected if 
the gap is created by dynamical friction. The differences for 
the subsamples with the smallest and largest gaps are ~ 1.3 
magnitude using the regular LF and ~ 2.8 magnitude using 
the relative LF. 

- The faint-end slope also shows a dependence on the mag¬ 
nitude gap, although the results are less significant. The 
Schechter faint-end slopes obtained for the four subsamples 
follow a trend moving from smaller to larger gaps. The a pa¬ 
rameter changes from a - -1.23 to a = -0.78 in the regular 
LF, and from a - -1.26 to a = -0.77 for the relative LF. We 
also fit an exponential to the data, finding the same trend in 
both cases. This is unexpected because the dwarf galaxy pop¬ 
ulation should not be affected by dynamical friction. Other 
processes, such as more radial orbits or early dwarf galaxy 
disruption, may play a role. 

The uncertainties in the faint-end slopes are mainly due to 
photometric uncertainties. To improve these results, it is impor¬ 
tant to study deep spectroscopic LFs for systems with large mag¬ 
nitude gaps. This is part of the future plans of the FOGO project. 
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